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TIME-DEPENDENT  SLOW-TIME-SCALE  THEORY 
OF  FREE-RUNNING  AND  PHASE-LOCKED 
GYROTRON  OSCILLATORS 

I .  INTRODUCTION 

There  is  currently  considerable  interest  in  the  development 
of  high  power  phase-locked  gyrotron  oscillators.  These  devices 
have  the  potential  to  combine  the  high  efficiency  and  power 
associated  with  oscillators  with  the  coherence  and  phase  control 
properties  associated  with  amplifiers.  Although  previous 
theoretical  work  on  gyrotron  operation[ 1-8 ]  has  been  successfully 
applied  to  the  development  of  cw  devices  for  heating  of  tokamak 
plasmas,  the  investigation  of  phase-locked  gyrotron  operation, 
particularly  involving  pulse-line  accelerator  driven  devices, 
presents  a  new  challenge  to  the  theorist.  These  devices,  based 
on  intense  relativistic  electron  beams,  are  characterized  by 
short  pulselengths  and  voltage  risetime  and  ripple  effects  which 
can  prevent  the  establishment  of  a  steady-state  equilibrium.  A 
consideration  of  non-stationary  operation  is  also  involved  in  the 
investigation  of  conditions  for  phase-locked  operation.  Time- 
dependent  effects  can  be  studied  using  a  particle-in-cell 
simulation  code  of  the  type  developed  by  Lin  and  co-wo r ke r s [ 9 ] . 

In  this  work  an  alternative  slow-time-scale  approach  is  used  to 
obtain  a  relatively  simple,  economical  analysis  of  the  time- 
dependent,  nonlinear  dynamics  of  pulsed  gyrotrons. 

The  time-dependent  theory  of  gyrotrons  has  been  considered 
by  Nusinovich  and  co-workers[ 10 , 11 ] ,  mainly  in  the  context  of 
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multimode  operation  and  mode  stability.  A  time  multimode  theory 
of  quasi-optical  gyrotrons  has  been  developed  by  Bondeson  et 
al.[12].  Early  work  on  mode  selection  and  phase-locking  of 
vacuum  tube  oscillators  was  carried  out  by  Van  der  Pol[13]  and 
Adler [14].  An  analytical  theory  of  the  conditions  for  phase¬ 
locking  gyrotrons  has  been  presented  by  Manheime r [ 1 5 , 16 ] .  The 
emphasis  in  the  present  work  is  on  the  modeling  of  transient 
effects  in  gyrotrons  whose  pulse  length  is  only  a  few  times 
longer  than  the  characteristic  risetime  of  cavity  fields.  The 
approach  is  developed  for  the  case  of  several  interacting  modes 
with  closely  spaced  frequencies  and  interacting  at  the  same 
harmonic  of  the  cyclotron  frequency.  Calculations  are  presented 
for  single  mode  operation  of  relativistic  gyrotron  with 
parameters  similar  to  an  experimental  device  under  test  at  the 
Naval  Research  Laboratory ( 17 ] .  Both  free-running  and  phase- 
locked  operation  are  considered.  The  phase-locking  signal  is 
introduced  via  a  prebunched  beam.  The  effect  of  a  voltage  ripple 
is  investigated.  The  possibility  of  phase-locked  operation  in 
the  presence  of  voltage  ripple  -  an  important  issue  for  proposed 
NRL  experiments  -  is  demonstrated.  Calculations  for  multimode 
operation  will  be  the  subject  of  future  work. 

The  time  theory  of  gyrotron  oscillators  derived  in  this  work 
is  based  on  slow-time-scale  equations  for  the  electron  motion 
similar  to  those  used  in  steady-state  models [ 7 , 18 , 19 ] .  Slow- 


time-scale  equations  for  the  cavity  rf  field  amplitude  are 
obtained  by  expressing  the  time-dependant  behavior  relative  to  a 
reference  frequency  w0  which  is  close  to  the  operating  frequency 
<t).  The  next  section  of  this  paper  contains  an  outline  of  the 
theoretical  approach.  The  third  section  contains  the  results  of 
calculations  and  the  last  section  presents  some  conclusions  drawn 


from  this  research. 


II.  THEORY 


Consider  a  gyrotron  with  a  cylindrical  resonator  and  a  thin 
annular  beam.  The  electrons  follow  helical  trajectories  in  the 
applied  axial  magnetic  field  about  guiding  centers  located  at  a 
radius  R0  from  the  symmetry  axis.  The  electron  beam  interacts 
with  one  or  more  resonator  modes,  which  are  assumed  to  be  TE 
modes  with  closely  spaced  oscillation  frequencies.  It  is 
convenient  to  look  at  time  effects  which  remain  after  a  reference 
frequency  w0  has  been  factored  out.  These  effects  are 
characterized  by  time  scales  which  are  much  longer  than  the  wave 
period  and  are  incorporated  in  time  mode  amplitudes  An ( t ) .  The 
total  transverse  electric  field  is  expressed  in  the  form 


Et  -  Z  An  ( t  )hn  (  z  )en  (  r  ,  9;  z  )exp< -ioo0  t )  (1) 

n  a  1 

where  en  is  a  transverse  mode  vector  function  defined  in  Ref.  18, 
and  hn ( z )  is  an  axial  profile  function  which  satisfies  the 
Helmholtz  equation; 


The  function  Pn»hnen  is  an  eigenfunction  of  the  empty  resonator 
in  the  single  transverse  mode  approximation  with  eigenf requency 
wno .  Thus  the  present  formulation  neglects  beam  loading  effects 
on  the  mode  axial  profile.  N  is  the  number  of  interacting  modes 
The  slowly  varying  mode  amplitudes  have  the  form: 

A„(t)  -  f„  <t)e-i*"  (t)  (3) 

The  phase  derivative  ^“d^/dt  corresponds  to  the  frequency 
difference  between  the  mode  oscillation  frequency  «n  and  «0 . 

The  transverse  electric  field  satisfies  the  wave  equation: 

.  32Et  3  Jt 

72  Et  _  A -  -  u0 -  (4) 

c2  3 t2  3 t 

where  Jt  is  the  transverse  ac  current  density,  c  is  the  speed  of 

light  and  fja  is  the  permeability  of  free  space.  MRS  units  are 
used  throughout  except  as  noted.  Substituting  Eq.(l)  into 
Eq.(4),  introducing  a  mode  Q  factor,  Qn ,  to  account  for 
diffraction  and  ohmic  losses,  and  noting  that  dAn/dt<<w0An  one 


obtains 


where  wno  is  the  empty  resonator  eigenf requency  for  the  nth  mode. 
Multiplying  Eq.(5)  on  both  sides  by  ( hn en ) * ,  integrating  over  the 
cavity  volume  and  applying  the  orthogonality  properties  of  the 
transverse  vector  functionsl 18  ]  leads  to  the  mode  excitation 
equation : 


r  d 

Wo 

-  ic2  fjo  _ 

+ - i(w0 

~  Wno  MAn  <  fc)  "  dzda 

Ldt 

2Qn 

J  2Nnco  i 

Nn 

is  the  axial 

mode  normalization  integral: 

iu>  t 


Nn  -  Pdz  |hn(z) 

J  Q 


To  evaluate  the  right  hand  side  of  Eq . ( 6 )  the  rf  current  density 
is  written  in  the  approximate  form(7]: 


J,.  =J  exp  (  -ioo  t ) 
4  “>0 


where 


J  =  f  d(wot)Jt  e1 
wo  J  n 


iw0  t 


Introducing  the  normalized  frequency  shift  6n  =  (  -oon  0  ) /t»>0  and 
time  t=<a)  t  and  substituting  Eq .  (  8  )  into  Eq .  (  6  )  leads  to: 


r  d  1  C  n 

—  + - i8n  An  (x)  - - dzda  h;  e;  •  J  e-  1  ’"n 

Ldx  2Q  J  2N„  m  .j 


To  calculate  the  ac  current  density,  the  interaction  with  the 
electron  beam  is  treated  in  the  single  particle  approximation.  A 
considerable  simplification  of  the  general  time-dependent  problem 
results  if  one  uses  the  fact  that  the  characteristic  rise-time  of 
fields  in  the  resonator  is  much  longer  than  the  electron  transit 
time  in  the  cavity.  In  this  case  one  can  use  a  quasi-steady- 
state  approximation  in  which  the  electron  trajectories  are 
calculated  for  rf  fields  with  fixed  amplitude,  .rj),  and 
linearized  phase,  =\|/n  '  (  x0  )  (  t-t0  )  .  The  slow-tiroe-scale 
nonlinear  electron  equations  of  motion  for  a  particular  harmonic 
are  readily  deduced  from  previous  steady-state  analyses [ 7 , 18 ]  and 
are  given  by: 


—  V"Js(knrt->  He 

ii  n>l  a 


k  +  i—  — 

L  rw.  dz 


_jne-i(  A+*n-(mn-s)2o  1 


(kn  rL 


u,  dh„ 


i  [  A+V>_ -(m  -s  )S„  ]  -  f, 

n  [  dhn 

z  Fn  Js  (kn  rL  )  Re  i— -  e 


where  ut«yvt/c  is  the  normalized  transverse  momentum  amplitude, 
u^-yv^/c  is  the  normalized  axial  momentum,  A  is  a  slowly  varying 
transverse  momentum  phase,  s  is  the  harmonic  number,  mn  is  the 
mode  azimuthal  index,  y  is  the  relativistic  mass  ratio,  kn  is  the 
mode  transverse  wavenumber,  rL  is  the  Larmor  radius  of  the  orbit, 
Js  (J's)  is  (the  derivative  of)  a  regular  Bessel  function,  5  is 
the  relativistic  cyclotron  frequency  correspondong  to  the  input 
beam  energy,  S0  is  the  orbit  guiding  center  azimuthal  angle,  and 
Fn  is  the  normalized  rf  field  amplitude: 


Fn  " 


J  (knRj 
mn  -s  n  ° 


m.c2  /[  n(  l-m^/x^2  )  ]  J  (  x^  ) 

mn 


(12) 


Quantities  with  a  "  have  been  normalized  according  to: 
i-z/ruo,  ?L-rL/rwo,  S-Crwo/c,  w0~«0rwo/c,  and  &n-knruo.  Ro 

denotes  the  orbit  guiding  center  radius,  e  is  the  electron 
charge,  m#  is  the  electron  mass,  x'n  is  a  zero  of  J'  for  the  mode 
n,  and  rwo  is  an  arbitrary  normalization  length  usually  chosen  to 
be  a  representative  cavity  wall  radius. 

The  ac  current  density  is  obtained  by  integrating  Eqs.(ll) 
for  an  appropriate  set  of  initial  conditions  at  the  cavity  input 
at  z0 .  For  a  cold,  phase-mixed  electron  beam:  ut(z0)=ut0, 
u2  (  z0  )  =»Uj,  0  ,  and  A(z0)=A0  is  uniformly  distributed  in  the  interval 
[ 0 , 2 n  ]  -  The  case  of  an  electron  beam  prebunched  by  passage 
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through  a  bunching  cavity  driven  at  a  frequency  wQ  can  be  modeled 
by  an  initial  phase  angle  distribution  of  the  form: 


A0  -  0O  +  qB  sin0o  -  *0  (13) 

where  0O  is  uniformly  distributed  in  the  interval  [0,2it],  qB  is 
the  bunching  parameter,  and  i0  is  an  initial  phase.  For  an 
annular  beam,  Eqs.(ll)  are  averaged  with  respect  to  S0 .  For  a 
thin  annular  beam  the  transverse  ac  current  density  is  given  by 


Substitution  of  Eqs.(14)  and  (9)  into  Eq.(10)  leads  to  the 
following  equations  for  the  mode  amplitude  and  phase: 


dr 


Re 


u 


*  -  -  i  [  A+4'  -(mn-s)S 


dzh;  < — J'  (  kn  rL  )e 


(15a) 


d*n 

dx 


6n  -  im 

I 

dzh;<—  J'(knrL)el[A+^  (m"  SUol>A  t 

v. 

uz  0  '  ”0  ) 

(15b) 


where  I  is  the  normalized  current  parameter: 


iMiMiwwum*unnuy 


mi  w  if  w  w 


M^c 


rwo  Jm-s(knRo} 


m0co>0  n(  1-m2  /x ' 2  ) J2 ( x '  )Nr 


(16) 


and  <  >a  denotes  an  average  with  respect  to  the  variable  a.  The 
time  behavior  of  a  free-running  or  driven  gyrotron  oscillator  is 
obtained  by  solving  Eqs.(lS)  together  with  (11)  subject  to 
appropriate  initial  conditions  as  discussed  above.  When  the 
phase  \|/n  of  a  driven  oscillator  (via  beam  prebunching)  converges 
to  a  constant  value,  the  oscillator  is  said  to  be  phase-locked. 
The  locking  frequency  is  &>0  in  the  present  formulation.  When  the 
interaction  occurs  with  only  a  single  circularly  polarized  mode, 
the  ac  current  density  averaged  with  respect  to  A0  is  independent 
of  the  guiding  center  azimuthal  position  so  that  the  average  with 
to  S  may  be  omitted. 
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III.  RESULTS 


Calculations  have  been  carried  for  a  high  voltage  35  GHz 
gyrotron  configuration  similar  to  the  NRL  experiment  recently 
reported  by  Gold  et  al.(17].  The  peak  voltage  and  current  of  the 
annular  beam  are  taken  to  be  650  kv  and  1.5  kA,  respectively. 

The  beam  guiding  center  radius  is  R0-1.16  cm,  the  cavity  radius 
is  1.6  cm,  and  the  longitudinal  profile  of  the  cavity  fields  is 
assumed  to  be  sinusoidal  with  length  4  cm.  The  operating  mode  is 
the  TEg  2  x  circular  mode  with  polarization  counter-rotating  to  the 
beam  rotation.  The  beam  pitch  ratio  asvt/v2=l.  The  cavity  Q 
factor  Q»250.  The  applied  magnetic  field  is  2.4  Tesla.  Spreads 
in  beam  guiding  center  and  pitch  ratio  are  neglected  as  are 
space-charge  effects.  In  these  calculations  the  cavity  is  seeded 
with  low  amplitude  rf  fields  (1  kV/cm  at  the  beam  in  the  present 
case).  The  results  of  calculations  for  free-running  and  phase- 
locked  oscillator  operation  are  presented  below.  The  effect  of  a 
voltage  ripple  is  shown. 

The  first  set  of  calculations  correspond  to  an  idealized  50 
nsec  voltage  pulse  with  a  3  nsec  rise  from  75%  of  peak  to  peak 
voltage  and  no  ripple.  The  beam  current  and  pitch  ratio  are 
assumed  to  vary  linearly  with  voltage  and  as  V1-33,  respectively, 
during  the  voltage  rise  up  to  the  maximum  values  given  above. 
Results  for  the  free-running  oscillator  are  shown  in  Figures  1 
(a)— 1(c)  .  These  Figures  show,  respectively,  the  electronic 


efficiency,  operating  frequency,  and  rf  phase  of  the  oscillator 
plotted  as  a  function  of  time.  Figure  1(b)  plots  the  shift  in 
operating  frequency  w  relative  to  <*>0  which  is  taken  to  be  the  the 
cold  cavity  eigenf requency  (which  is  close  to  35  GHz).  These 
plots  indicate  that  steady-state  operation  is  achieved  after 
about  8  nsec.  The  steady-state  efficiency  is  17.7%, 
corresponding  to  an  output  power  of  173  MW,  and  the  steady-state 
frequency  shift  5r  R «( ajr  R -co0  )/w0  due  to  beam  loading  is  5.5xl0-4 
which  corresponds  to  19.2  MHz.  This  frequency  shift  leads  to  a 
linearly  increasing  rf  phase  with  time  as  shown  in  Figure  1(c). 

Results  for  an  oscillator  driven  via  a  prebunched  beam  are 
given  in  Figures  2(a)-2(c).  The  beam  is  modulated  at  a  frequency 
oo0  and  the  modulation  amplitude  is  modeled  by  the  bunching 
parameter  qB=0.5.  An  estimate  for  the  locking  signal  power 
required  to  obtain  a  beam  with  qB»0.5  is  given  below.  In  this 
calculation  the  difference  between  the  normalized  locking  signal 
frequency  and  the  normalized  free-running  oscillator  (FRO) 
frequency  is  SL  =  (  w0  -<*>F  R  )  /oo0  =1 . 0  5  xlO~  3  or  Af=35  MHz.  Figures  2(b) 
and  2(c)  show  that  phase-locked  operation  is  obtained  after  30 
nsec  into  the  voltage  pulse.  The  phase-locked  condition  is 
indicated  by  a  frequency  shift  oo-wo=0  and  constant  rf  phase.  As 
shown  in  Figure  2(a),  in  this  case  the  oscillator  efficiency  is 


reduced  slightly  to  16%  while  phase-locked. 


If  the  difference  between  the  locking  signal  frequency  and 
the  FRO  frequency  Af  is  too  large,  phase-locking  will  not  occur. 
Instead  a  periodic  perturbation  of  the  oscillator  behavior  is 


observed  as  shown  in  Figures  3(a)-3(c).  The  locking  signal 
detuning  in  this  case  is  SL»1.55xl0-3  or  Af-54  MHz.  The  beat 
frequency  of  the  perturbation  is  about  28  MHz.  The  non-phase- 
locked  behavior  of  the  rf  phase  is  clearly  shown  by  Fig.  3(c). 

The  case  of  a  FRO  with  a  voltage  ripple  is  shown  in  Figures 
4(a)-4(c).  The  voltage  ripple  amplitude  is  2%  and  the  ripple 
period  is  10  nsec.  Other  parameters  are  the  same  as  for  the  FRO 
case  discussed  above.  Figures  4(a)  and  4(b)  show  that  this 
voltage  ripple  causes  a  3.4%  amplitude  ripple  in  the  efficiency 
and  a  0.06%  amplitude  ripple  (21  MHz)  in  the  frequency.  Figures 
5(a)-5(b)  show  the  results  of  driving  this  oscillator 
configuration  with  a  prebunched  beam  with  q„=0.5  and  modulation 
frequency  shift  &L=5.5xlO_4  (19  MHz).  Figures  5(b)  and  5(c)  show 

that  phase-locked  operation  is  achieved  in  the  sense  that  after 
about  30  nsec  the  oscillation  frequency  averaged  over  the  ripple 
period  equals  the  drive  frequency  and  the  averaged  rf  phase  angle 
is  a  constant.  As  in  the  FRO,  the  voltage  ripple  leads  to  a 
0.06%  frequency  ripple  in  the  PLO.  The  frequency  ripple  causes  a 
corresponding  13°  amplitude  ripple  in  the  phase  angle.  Figures 
6(a)-6(c)  represent  the  case  of  an  unphase-locked,  driven 
oscillator  with  a  2%  voltage  ripple.  The  drive  signal  frequency 
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shift  is  Sl-2.5x10“3  (88  MHz).  The  beat  frequency  for  the 
perturbed  oscillator  is  about  34  MHz.  The  failure  to  achieve 
phase-locked  operation  is  clearly  shown  by  the  phase  angle 
plotted  in  Figure  6(c). 

To  estimate  the  input  power  needed  to  obtain  a  beam  bunching 
parameter  qB«0.5  for  circularly  polarized  radiation  and  a 
gaussian  rf  profile  in  the  bunching  cavity,  one  can  apply  the 
small  signal  gyroklystron  theory  developed  by  Tran  et  al  (20]. 
The  required  input  power  in  Watts  is: 


Pin(W) 


'I  nc3  m2  (x^2-m2)  J£  (x'n  )  y|6?03Z0  q|  e‘ 
2 ’I  2//0  e2  J2  (knR0)  nx  (4  3  Q1 


(17) 


where  /y1-n02to  L1/(3zoX  is  the  normalized  interaction  length  of  the 
bunching  cavity,  uA  is  the  normalized  drift  length,  Qx  is  the 
bunching  cavity  Q  factor,  and  x= ( wQ -B ) Lx /4 vz o  the  bunching  cavity 
magnetic  field  detuning  parameter.  The  other  parameters  are 
defined  above.  Typical  choices  for  /jx  and  /ud  are  2  and  3, 
respectively.  Assuming  the  same  mode,  magnetic  field  and  beam 
parameters  in  the  bunching  and  power  cavities,  and  Q1=200,  one 
obtains  Pin=«530  kw  as  the  power  required  to  phase-lock  a  gyrotron 
oscillator  with  an  output  power  of  over  160  MW  for  an  effective 
gain  of  over  24  dB . 


IV.  CONCLUSIONS 


A  time  slow-time-scale  theory  and  a  computer  code  for 
numerical  calculations  have  been  developed  to  investigate  pulsed 
free-running  and  phase-locked  gyrotrons.  This  code  is  expected 
to  be  a  valuable  resource  for  the  design  and  interpretation  of 
planned  phase-locked  gyrotron  experiments.  The  calculations 
presented  here  show  the  feasibility  of  observing  phase-locked 
operation  with  50  nsec  pulses.  Phase-locking  over  a  significant 
bandwidth  is  shown  to  be  possible  with  a  locking  power  24  dB  less 
than  the  oscillator  output  power  for  a  single  prebunching  cavity. 
The  effect  of  voltage  ripple  has  been  investigated  and  shown  not 
to  prevent  phase-locked  operation,  however,  applications 
requirements  for  rf  phase  control  may  place  significant  limits  on 
allowable  voltage  ripple. 
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